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ABSTRACT 

Aims. In this work the latitude dependent stellar spot rotation is investigated based on dynamo models, and using similar analysis 
techniques as in the case of real observations. The resulting surface differential rotation patterns are compared to the known input 
rotation law used in the calculation of the dynamo model. 

Methods. Snapshots of the dynamo simulations are used to measure the surface differential rotation. The maps of the magnetic 
pressure at the surface are treated similarly to the temperature maps obtained using Doppler imaging techniques, and a series of 
snapshots from the dynamo models are cross-correlated to obtain the shift of the magnetic patterns at each latitude and time point. 
These shifts are used to study the surface rotation pattern over a wide latitude range at different epochs during the activity cycle, and 
are compared to the known input rotation law. 

Results. Two different rotation laws are investigated, one solar-type law and one with axis distance dependent rotation. Three different 
dynamo calculations are carried out based on the axis distance dependent law: one with only large scale dynamo field, one with 
additional strong small scale field and one with weaker small scale field. The surface differential rotation patterns obtained from the 
snapshots of all the four dynamo calculations show variability over the activity cycle. Clear evolution and variation in the measured 
surface rotation patterns is seen, but in the models using only the large scale dynamo field the measured rotation patterns are only at 
times similar to the input rotation law. This is due to the spot motion being mainly determined by the geometric properties of the large 
scale dynamo field. In the models with additional small scale magnetic field the surface differential rotation measured from the model 
follows well the input rotation law. 

Conclusions. Here, for the first time, the surface differential rotation patterns are investigated in detail based on dynamo calculations. 
The results imply that the stellar spots caused by the large scale dynamo field are not necessarily tracing the stellar differential rotation, 
whereas the spots formed from small scale fields trace well the surface flow patterns. It can be questioned whether the large spots 
observed in active stars could be caused by small scale fields. Therefore, it is not clear that the true stellar surface rotation can be 
recovered using measurements of large starspots, which are currently the only ones that can be observed. 

Key words, magnetic fields - MHD - Stars: activity - Stars: rotation - starspots 

1. Introduction tion zones predict solar-type DR, where the equator is rotating 

faster than the poles. However, Kitchatinov & Riidiger (2004) 
Differential rotation (DR) is one of the key ingredients in the have shown that anti-solar DR, which wa s sugg ested by obser- 
stellar dynamo models. Together with the helical turbulence and vations of several active stars (e.g., Weber |2007 1, could arise as 
meridional flow it is responsible for the main features of t he so- a result of intensive meridional circulation, 
lar and stellar magnetic activity (see, e.g., Riidiger et al. |1986l 
b ■ Choudhuri et al. [1995] Brun & Toomre 2002 1. The surface DR Stellar differential rotation is usually characterised by the 
- -' of the Sun has been known for a long time from observing the surface shear 5Q - Q eq - ii pole , where Q eq and Q pole are the 
rotation of sunspots and other magnetic elements at different lat- rotatlon rates at the ec l uator and the P oles ' respectively. The sur- 
itudes (e.g., Balthasar et al. 1986 ). The rotation period at the face shear 18 related to the la PP ln 8 tlme ' ^ = 2n / m : For stars ' 
solar equator is approximately 30% shorter than the period at a surface rotation law of the form Q = Q eq (l - fcsin 0) is usu- 
the poles. Helioseismology has revealed that this pattern per- all Y assumed, where 6 is the latitude. The measurements of sur- 
sists t hroughout the ent ire convection zone (e.g., Thompson et al. face P R on r apidly rotating cool stars is summarised by Barnes 
IT9961 Coo et al. [19981 ), while the radiative core rotates rigidly, et al. (|2005>. These measurements seem to show an increase 
Between the core and the convection zone there is a transition in the magnitude of differential rotation towards earlier spec- 
layer^ tachocline, with a strong radial shear (Spiegel & Zahn tral tyP es ' which is consi stent with the theoretical calculations of 
11992) . Riidiger & Kiiker (2002). Recently, a comparison of the models 

™ , , , . ... , . . j j ( and observations was carried out by Kiiker & Riidiger (2008 1. 
Theoretical calculations predict that across the Hertzsprung- J ° * ' 

Russel diagram DR becomes stronger with increasing rotation Intriguingly, temporal evolution of the surface DR has been 

period and evolutionary status (Kitchatinov & Riidiger 1999). reported for two young single K stars, AB Dor and LQ Hya (see, 

In general, models for global circulation in outer stellar convec- e.g., Donati et al. 120031 Jeffers et al. 12007) . Donati et al. (2003) 



2 



H. Korhonen and D. Elstner: Investigating variation of latitudinal stellar spot rotation 



hypothesise that these temporal variations could be caused by 
the stellar magnetic cycle converting periodically kinetic energy 
within the convective zone into large-scale magnetic fields and 
vice versa, as originally proposed by Applegate (11992) . 

Measuring stellar differential rotation is not straight forward. 
Stars are faint and apparently small, thus only the closest su- 
per giants can be directly spatially resolved with the current in- 
struments and methods (e.g., Gilliland & Dupree 1996; Monnier 
et al. 120071) . Currently the best way to study stellar surface DR 
in detail, both the strength and the sign, is by using surface 
maps obtained with Doppler imaging technique (see, e.g., Vogt 
et al. |1987l Piskunov et al. |19901 ). One way, similar to using 
sunspots and other magnetic features to study the solar DR, is to 
use several of the surface maps and cross-correlate them to inves- 
tigate the changes in the locations of the spots (see, e.g., Barnes 
et al. [20001 Weber et al. 120531 Kovari et al. |20Q7T >. The other 
much used way is so-called ^-landscape technique in which 
a DR law is implemented into the Doppler imaging code, and 
many maps are obtained to study which parameters give the best 
soluti on (e. g., Petit et al. 120021 Marsden et al. [2006 ; Dunstone 
et al. 120081 1. Stellar surface DR has also been studied by com- 
bining spot latitude information from Doppler images and spot 
rotation period from contemporaneous photometry (Korhonen et 
al. 2007) and by using Fourier transforms to study shapes of the 
spectral line profiles (see, e.g., Reiners & Schmitt l2003) . 

All the methods for measuring stellar surface DR using 
Doppler imaging results suffer from the often restricted latitude 
range the starspots occur on, and from possible artifacts in the 
maps due to systematic errors. These systematics can for exam- 
ple rise from incorrect modelling of the spectral line profiles. 
Additionally, in ^-landscape technique a predefined rotation 
law is assumed. On the other hand, in cross-correlation the time 
difference between the maps is crucial, with too small difference 
no change has time to occur, and with too long difference spot 
evolution due to flux emergence and disappearance could have 
occurred. Therefore, it is very demanding to obtain reliable mea- 
surements of stellar surface differential rotation, and currently it 
is basically impossible to obtain information on the internal ro- 
tation of stars other than the Sun. 

In this study we investigate the surface differential rota- 
tion obtained from snapshots of dynamo models using cross- 
correlation methods. These snapshots provide a very good tem- 
poral base for studying the measured surface differential rotation 
at different times over the stellar activity cycle. The obtained sur- 
face DR values are compared to the known internal rotation law 
used in the dynamo modelling. 



2. Dynamo model and analysis methods 

2.1. Model 

The model used here was first developed to theoretically study 
the flip-flop dynamo behaviour. In flip-flop phenomenon the 
spots concentrate on two permanent active longitudes that are 
separated by 180°, and the main spot activity switches be- 
tween the longitudes every few years (e.g., Jetsu et al. 119931 
Berdyugina & Tuominen 1998; Jarvinen et al. 120051 . 

The dynamo is modeled with a turbulent fluid in a spheri- 
cal shell. Two different rotation laws are chosen for investiga- 
tion: one similar to the solar one but with a smaller difference 
between core and surface rotation (LI), and another one with 
axis distance dependent rotation (L2), similar to the models pre- 
sented by Moss (120041 ). The mean electromotive force contains 
an anisotropic alpha-effect and a turbulent diffusivity. The non- 
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Fig. 1. Illustration of the internal rotation laws used in this study, 
(a) The solar-type rotation law with the equator rotating faster 
than the poles (LI), (b) Axis distance dependent law, where parts 
furthest away from the axis are rotating the fastest (L2). 

Table 1. Dynamo models. The name of the model, radius of the 
inner boundary in stellar radii, dynamo numbers, the distribution 
of the magnetic energy in different modes, and the cycle period 
are given. 



name 


fin 


c n 


r 


eo 


e\ 


e„ 


period 




[ R*] 






[B cq ] 


[Be,] 


[Be,] 


[years] 


LI 


0.4 


184 


20 


0.06 


0.05 


6 ■ 10-' 


6.8 


L2 


0.27 


250 


13 


0.12 


0.16 


0.008 


2.7 


L2a 


0.27 


250 


13 


0.10 


0.2 


0.04 


2.7 


L2b 


0.27 


250 


13 


0.11 


0.15 


0.01 


2.7 



linear feed-back of the magnetic field acts on the turbulence only. 
The boundary conditions describe a perfect conducting fluid at 
the bottom of the convection zone and at the stellar surface the 
magnetic field matches the vacuum field. The model is discussed 
in detail by Elstner & Korhonen (2005 ) and Korhonen & Elstner 
d2005l) . 

The rotation laws used in this study are illustrated in Fig. [TJ 
In LI the internal rotation is solar-type, i.e., the equator is rotat- 
ing faster than the polar regions. The inner boundary of the con- 
vective zone is at the radius 0.4 R*. L2 is an axis distance depen- 
dent rotation law, where the parts furthest away from the rotation 
axis are rotating the fastest. The areas close to the axis, i.e., on 
the surface the polar regions, do not rotate at all. In this model the 
inner boundary of the convective zone is at the radius 0.27 R*. 
The dynamo numbers C a = a * R^/tjt and Cq = AQ * R+/t]t, 
and the distribution of magnetic energy (eo in the axisymmetric 
mode, e \ in the mode with azimuthal wavenumber m = 1 , and e„ 
in the higher modes) are given in Table Q] 

Both models lead to a mixed field of a non-axisymmetric 
and an oscillating axisymmetric mode. The non-axisymmetric 
modes rotate rigidly with Q m = -0.32 degrees /day for LI and 
Cl m = -2.9 degrees/day for L2. It corresponds to the surface 
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Fig. 2. Illustration of the clockwise spot motion (circles) during 
the flip-flop event. Shown is the surface magnetic pressure of 
the non-axisymmetric mode. The circles along latitude mark the 
maximum of the poleward migrating axisymmetric field at dif- 
ferent times. Spots appear only at one spiral, where the signs of 
the two modes are equal. Note, that the stellar rotation and also 
the magnetic pattern are anti-clockwise. 



rotation at about 55° latitude for LI and 65° for L2. The oscil- 
lation period for the full cycle of the axisymmetric mode is 6.8 
years for LI and 2.7 years for L2. In Fig. |2] the surface mag- 
netic pressure of the non-axisymmetric mode for L2 is shown. 
The schematic rings demonstrate the upward motion of the ax- 
isymmetric field belt of one polarity during half of a cycle. The 
small circles indicate the spot motion due to the field geome- 
try (in the coro taring frame of the field pattern). The spots move 
during one third of a half cycle (40° to 70°) about 100° opposite 
to the stellar rotation. Therefore we expect nearly a cancella- 
tion of the spot rotation during the flip-flop cycle. The detailed 
behaviour depends on the spiral pattern of the field and the mi- 
gration speed of the axisymmetric field. A similar explanation 
for the differential rotation of active longitudes on the sun was 
already discussed by Berdyugina et al. (2006 1. 

In order to study the effect of small scale spots additional 
simulations were run for the model L2, where small poloidal 
field loops were randomly added into the convection zone once 
per day (laws L2a and L2b). The magnetic energy due to the 
additional field injection in model L2a is about half of the en- 
ergy of the axisymmetric field component in the quasi-stationary 
regime. The fields appear as localised spots on the surface, which 
move with the surface rotation. With that strong field injection 
(L2a) the field structure of the global dynamo mode is nearly 
hidden by the small scale spots. Therefore in a second run 10 
times weaker field loops (model L2b) were added, which led to 
a final energy of the small scale component of only l/100th of 



the axisymmetric field. Here the global field again dominates the 
surface spot structure. One could also correlate the small scale 
field injection with the large scale magnetic field, but for sim- 
plicity a totally uncorrelated field injection was chosen. 

2.2. Measuring the surface differential rotation 

When using the dynamo models to study the surface differen- 
tial rotation, the temporal resolution can be chosen high enough 
that standard cross-correlation methods can be used. Snapshot 
maps of the surface magnetic pressure are taken from the dy- 
namo calculations and treated as representations of the magnetic 
structures at that time point. These maps can be treated the same 
way as the temperature maps obtained from Doppler imaging, 
therefore the same techniques as in the case of the real obser- 
vations can be used to analyse the model snapshots. Magnetic 
pressure maps have been taken from 36 time points over the ac- 
tivity cycle. These time points are relatively close to each other: 
separated by about 50 days in the LI, and by about 18 days with 
the L2. Examples of snapshots for LI and L2 are shown in on- 
line Figs.[3]&S respectively. 

In the analysis each time point is cross-correlated with the 
following one and the behaviour is studied from equator to the 
visible pole. The dynamo model has a grid of 43 x 43 surface 
points. This makes the longitudinal resolution of 8.4° and latitu- 
dinal of 4.2°. The equator is going through in the middle of the 
22nd latitude strip. The cross-correlations are obtained for each 
latitude between the latitudes centred at 4.2° and 83.7°. The last 
latitude strip closest to the pole is not used due to the lack of sig- 
nal. The measured shift in degrees/day at each latitude is com- 
pared to the expected shift obtained from the rotation law used 
in the dynamo model. The spots are migrating in the dynamo 
simulation, to remove this field migration all the measured shifts 
are normalised to the shift at the lowest latitude used in the in- 
vestigation. 

3. Results 

Examples of results from cross-correlating the 36 maps for L 1 
are shown in Fig. [5] The plots give the shift in degrees/day 
for each latitude on the visible hemisphere (crosses). The last 
plot shows the average of the measurements from all the cross - 
correlations, with standard deviation of the measurements as the 
error. In the plots the dashed line is the input rotation at the stel- 
lar surface. 

The behaviour of the measured spot rotation clearly changes 
from cross-correlation to cross-correlation, indicating that the 
appearing surface rotation pattern changes during the activity cy- 
cle. Furthermore, the measured surface rotation can take numer- 
ous different forms: typical solar-type rotation law (e.g., cross- 
correlation 11), almost no difference between the latitudes (e.g., 
cross-correlation 01) and even completely anti-solar rotation law 
with the pole rotating faster (e.g., cross-correlation 25). In gen- 
eral it is evident that in the case of LI the measured surface dif- 
ferential rotation at the low latitudes often is what one would 
expect from the models, but at higher latitudes the correlation is 
very poor. 

The examples of cross-correlation results from L2, shown in 
Fig. [6] also exhibit large variations in the surface rotation pat- 
terns over the activity cycle. In some cases the measured surface 
rotation at the lower latitudes is similar to the input rotation (e.g., 
cross-correlation 14), at times the surface rotation at the equator 
and at the polar region is similar to the input rotation, but the 
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Fig. 5. Examples of the cross-correlation results using the 36 individual snapshots from the dynamo model LI. The crosses give 
the shift in degrees/day for each latitude between 4.2° and 83.7°. The dashed line in the plots shows the input rotation law for the 
dynamo model at the surface. The last plot is the average behaviour obtained from all the snapshots with using the standard deviation 
of the measurements as the error. 



exact shape of the rotation curve is different from the expected 
(e.g., cross-correlation 16), but in most cases the rotation shows 
much less latitudinal variation than expected based on the input 
rotation. 



4. Discussion 

4.1. Surface differential rotation patterns over the activity 
cycle 

4.1.1. Model L1 

The surface DR measurements for model LI, shown in Fig. [5] 
cover slightly more than one activity cycle, which also includes 
a flip-flop event. Similar phase in the cycle as in the map 1 is 
reached in the map 28. Clear evolution in the surface rotation 
patterns are seen throughout the cycle. 

In general the shape of the surface rotation pattern over the 
latitude can be divided into four parts. Closest to the equator the 
surface rotation follows the internal rotation law, and in most 
maps three abrupt changes in the surface rotation occur at higher 
latitudes. Clear evolution in the latitude at which the breaks oc- 
cur is seen over the activity cycle. When the cycle advances the 
breaks drift towards the higher latitudes. Even at the point that 
the spots move to lower latitude, the breaks continue drifting 
towards the polar regions. The highest latitudes for the breaks 
are obtained just before the time when the second spot starts to 
appear in the flip-flop event. When the breaks abruptly move to- 
wards the equator again, the rotation measurements for the lower 
latitudes become less ordered and appear chaotic. 

In this model the axisymmetric dynamo wave is not only a 
poleward migration as indicated in Fig. [2] It starts at the bot- 
tom of the convection zone at 30° latitude and expands out- 
wards in radius and latitude. For lower latitudes the expansion 
along the non-axisymmetric spiral gives now a counter clock- 
wise acceleration. This leads to the qualitative agreement with 
the stellar rotation law near the equator. The effect depends on 
the expansion speed and the pitch angle of the underlying non- 



axisymmetric field. For the rotation law dependent small pitch 
angle near the equator this is a huge effect. At higher latitudes the 
spot motion is determined by the field geometry and the migra- 
tion speed of the poleward moving dynamo wave and therefore 
the rotation is reduced by the clockwise motion of the spot along 
the spiral. At the very polar region the wave moves slower and 
the rotation appears again faster. This behaviour depends on the 
phase of the cycle. If the field reversal between two cycles lies at 
mid latitudes (50°) we see the rotation of the non-axisymmetric 
field, which coincides here with the stellar rotation (cf. Cross- 
correlation 1 to 5 and 26 to 30). 

4.1.2. Model L2 

The measurements of the surface DR using model L2 shown in 
Fig.|6]cover almost two activity cycles. The same epoch as in the 
first cross-correlation is again obtained around the time of the 
cross-correlation 18. Evolution in the surface rotation patterns 
are seen, but they are not as clear as in the case of LI . 

For most epochs the surface rotation is similar at all the 
latitudes, showing relatively small differential rotation. Often 
two areas with distinctly different rotation behaviour are seen, 
though. One where the magnetic surface patterns lag behind, and 
one where they rotate faster than the general surface rotation. 
The mainly constant spot velocity in this model goes back onto 
much more clear poleward migration of the oscillating mode be- 
cause of the different shape of the inner rotation law. Here we 
are missing the equatorward motion of the axisymmetric mode 
and do not see any correspondence to the input rotation law. The 
small positive values sometimes near the equator may be some 
similar, but much weaker, effect as for the model LI. The neg- 
ative jumps migrating poleward are connected to the rigid ro- 
tation of the non-axisymmetric dynamo mode, seen where the 
poleward migrating axisymmetric mode is about zero. 

The DR measurements for L2 usually show very rigid ro- 
tation. During the activity cycle the most rigid rotation occurs 
around the time of the flip-flop event, and the best fit to the rota- 
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Fig. 6. The same as Fig. [5] but using L2. Note that occasionally very large negative shifts are measured and one point can fall outside 
the plotting range. 



tion law around the equator and the polar region is obtained just 
before the flip-flop event starts. Still, even at these time points 
the shape of the measured rotation law is very different from the 
shape of the input rotation law. Two deviations from the gener- 
ally rigid rotation are often seen. In many cases a small latitude 
range shows larger than normal positive shifts. This area is re- 
stricted to latitudes >20°. Also, an area of negative shifts is seen. 
This area migrates polewards with the general migration of the 
magnetic elements. At the time leading to a flip-flop event this 
area occurs in between the two latitudinal spot belts. 



4.2. Small scale spots versus large scale dynamo spots 

It is interesting to note that the DR measured from the models 
only including large scale dynamo field follows the input DR 
the best at the latitudes close to the equator, and breaks down 
completely closer to the poles where the magnetic pressure is 
the highest. This could imply that the strongly magnetic regions 
do not follow the differential rotation as well as the less magnetic 
areas. 

On the first glance a correspondence of the large scale dy- 
namo field with the local differential rotation cannot be expected. 
Only the non-axisymmetric features lead to an observable rota- 
tion. The main component here is a large scale field with az- 
imuthal wavenumber m — 1. This field rotates rigidly with 
a value between the minimum and maximum rotation speed. 
Because of the faster rotation at the equator this field is shaped 
like a leading spiral. Together with the oscillating poleward mi- 
grating axisymmetric field mode a counter clockwise spot mo- 
tion is introduced, which may be latitude dependent (cf. Fig. [2]). 
All these effects have no direct connection to the differential ro- 
tation of the stellar surface. Therefore, the rotation law cannot be 
recovered by measuring the spot motion of the large scale field. 
On the other hand, at least for LI, a good agreement of the cross 
correlation result with the actual rotation law is found at low lati- 
tudes. This is a consequence of the equatorward extension of the 
field. Here, a connection to the stellar rotation law is obtained 
through the pitch angle. 



Tests were carried out whether the DR patterns can be re- 
covered from models including also small scale fields, by adding 
these fields to the model L2. Two different strengths for the small 
scale fields were used: strong field (L2a) and a 10 times weaker 
field (L2b). For more details on the models see Section |2~T1 and 
TableQ] 

The examples of snapshots for L2a and L2b are shown in on- 
line Figs. [7]&[8j respectively. In the analysis 99 maps are used. 
In these models, especially in L2a, which has strong small scale 
field, the changes are so rapid that the time step between the 
snapshots is now reduced to approximately 6 days. This means 
that with the 99 snapshots we cover more or less the same time 
period as with the 36 snapshots earlier. 

Examples of the cross-correlation results from L2a are 
shown in Fig. [9] The cross-correlation is done for 99 snapshots of 
magnetic pressure, but only every other cross-correlation curve 
for the time period covering a bit more than one flip-flop cycle is 
shown in the plot. Here, the latitude 12.6° is used as the reference 
latitude due to the anomalous shifts occasionally measured at the 
lower latitudes. As can be seen, now the recovered surface DR 
follows well the input rotation. Still, some time points of anoma- 
lous rotation are seen, but, as shown by the average rotation, in 
general the rotation is well reproduced. 

The results for L2b are shown in Fig. [10] Now the latitude 
4.2° is again used as the reference latitude. With a weaker small 
scale field of the model L2b the results are between the results 
from the original L2 and from L2a, in which stronger small scale 
field was used. Again, the general rotation law is recovered, but 
the high latitudes start to again show almost rigid rotation. 

Addition of the small scale fields to the dynamo model en- 
ables us to recover the input rotation law. But it remains an 
open question if the observable large spots can be attributed to 
small scale magnetic features. Recently, Lehtinen et al. (1201 II) 
found active longitudes rotating significantly slower than the 
star, which according to the current investigation, favors spots 
as signs for the large scale field. 
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Fig. 9. Examples of the cross-correlation results using all the 99 individual snapshots from the dynamo model L2a. The symbols are 
the same as in Fig. [5] The last plot is the average behaviour obtained from all the 99 snapshots with using the standard deviation of 
the measurements as the error. Latitude 12.6° is used as the reference latitude. 
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Fig. 10. As for the Fig. [9] but now for L2b and using latitude 4.2° as the reference latitude. 



4.3. Comparison to observations 

The results presented here clearly show that the rotation law ob- 
tained from dynamo models with large scale fields breaks at the 
latitudes where most of the spots occur. 

A literature search was carried out to investigate whether 
or not this phenomenon is also seen when cross-correlating 
Doppler images. One significant difference between the model 
calculations used here and the real Doppler images, is that 
when cross-correlating the models there is always signal, even 
if the signal is weak, i.e., not resulting in a starspot, whereas in 
Doppler images the latitudes devoid of spots only result in noise 
in the cross-correlation. 

The literature search reveals that at times the cross- 
correlation result reproduces very well the sin 2 -type surface ro- 
tation law within observational errors, which at times can be 



very large. A relatively good fit is seen for example in RS CVn- 
type binarie s IM P eg (Weber et al. 120051 1 and IL Hya (Weber & 
Strassmeier 11998) . young single K dwarfs AB Dor (Donati & 
Collier Cameron \TWj\ and LQ Hya (Kovari et al. 120041, and 
for a post T Tauri star RX J1508.6 + 4423 (Donati et al. l2000l) . 
On occasions the sin 2 -law has the same general shape as what 
is obtained from cross-correlating latitude strips in the Doppler 
images. For example the shifts measured for a late-type sin- 
gle star PZ Tel follow well the general shape of the sin 2 -law, 
but there are latitudes, especially around the equator, where the 
measurements deviate significantly from the sin 2 -law (Barnes et 
al. 12000). Also for VI 192 Ori the general shape of the rotation 
law is reproduced, but deviations are seen at the higher latitudes 
(Strassmeier et al. 2003). On KU Peg the latitudinal rotation 
rates form a complex double peaked pattern, where the equator 
and the latitudes around 45° rotate faster than the general rota- 
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tion (Weber & Strassmeier IIM] Weber et al. 120051 . There are 
also examples of targets where Doppler images obtained from 
different spectral lines give slightly different latitudinal rotation 
behaviour, e.g., <x Gem (Kovari et al. l20071 > and UZ Lib (Vida et 
al. 120071 

One has to also note that the results on cross-correlating 
Doppler images do not commonly extend to higher latitudes 
than 60°. Often this is due to the polar spot, which is seen in 
many of the maps. This constant feature would not result in shifts 
in spot configuration between two consecutive Doppler images, 
and thus it would not give any signal in the cross-correlation. 
Still, this feature is often large, and the strongest spot on the sur- 
face, i.e., analogous to the maximum magnetic pressure in the 
dynamo models, and exactly the location where the rotation be- 
haviour is seen to deviate from the input rotation. Thus, it is very 
difficult to conclude based on the Doppler images whether or not 
the observations also show the behaviour where the surface ro- 
tation in strongly magnetic features deviates from the general 
rotation. 



5. Conclusions 

Starspots are generally used as tracers of the stellar surface dif- 
ferential rotation. If starspots are direct signs of the global stel- 
lar magnetic field this assumption does not seem to be valid, 
because in this case the spot motion would mainly be deter- 
mined by the geometric properties of the large scale dynamo 
field. Based on the dynamo calculations presented here, the sur- 
face differential rotation can only be recovered by following the 
spot motion of small scale spots frozen into the flow. 

In the dynamo models presented here the usual assump- 
tion of a positive alpha in the northern hemisphere is made, 
implying a poleward migration of the axisymmetric dynamo 
wave. This gives together with the leading spiral of the non- 
axisymmetric field geometry a strongly reduced or even counter 
rotating spot motion. This could possibly explain the reported 
detections of anti-solar differential rotation in some stars (e.g., 
Vogt et al. 119991 Weber et al. 120051 . Nevertheless, for equator- 
ward drifting dynamo waves, as in the Sun, the spot motion will 
more or less coincide with the surface differential rotation. In 
our model LI it coincides even quantitatively in the equatorial 
region. This may be due to a connection of the strength of the 
differential rotation to the pitch angle of the non-axisymmetric 
mode. The fact that mainly solar type rotation laws have been 
observed (see, e.g., Barnes et al. 20051 would imply that the 
dynamos in general show an equatorward drift. Further mod- 
els with more solar type dynamos are needed. In case that the 
observed spots are more local phenomena we have a direct cor- 
respondence of the spot motion to the stellar rotation law. This 
was nicely reproduced using the cross-correlation method. 

Here, stellar surface differential rotation was for the first time 
investigated in detail based on dynamo calculations. The follow- 
ing conclusions can be drawn from this study: 

1. Cross-correlation is a powerful method to follow differen- 
tial rotation using starspots as tracers, if the time difference 
between the maps being cross-correlated is appropriate for 
detecting the differential rotation induced changes. 

2. Spots due to the large scale dynamo field are not necessarily 
tracing the real surface differential rotation. 

3. Equatorward drift of the dynamo field gives at least qualita- 
tively an agreement with the differential rotation, as was seen 
for the model with a solar type rotation law (LI). 



4. The differential rotation obtained from the spot motion varies 
during the dynamo cycle. 

5. Spots caused by small scale magnetic fields are good tracers 
of the stellar differential rotation. 

6. Comparing direct tracers of differential rotation (e.g., 
analysing shapes of spectral lines) and obtaining the spot 
motion from cross-correlation can probably clarify the na- 
ture of the spots, i.e, whether they are caused by small scale 
field or are signs of the large scale dynamo field. 
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Fig. 3. Examples of snapshots from the dynamo model LI. In total 36 individual magnetic pressure maps covering a slightly more 
than one flip-flop event are used in the analysis. 
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Fig. 4. The same as for Fig. [3] except now for model L2. 
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Fig. 7. Examples of the snapshots from the dynamo model L2a. In total 99 magnetic pressure maps are used in the analysis. 




Fig. 8. The same as for Fig. [7] except now for model L2b. 



